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ABSTRACT 

We use the 2011 Rossi X-ray Timing Explorer (RXTE) proportional counter array 
(PCA) data of the 401 Hz accreting pulsar and burster IGR J17498-2921 to perform 
timing analysis and time-resolved spectroscopy of 12 thermonuclear X-ray bursts. We 
confirm previously reported burst oscillations from this source with a much higher 
significance (8.8cr). We notice that the bursts can be divided into three groups: big 
photospheric radius expansion (PRE) bursts are about ten times more luminous than 
medium bursts, while the latter are about ten times more luminous than small bursts. 
The PCA field-of-view of these observations contains several known bursters, and 
hence some of the observed bursts might not be from IGR J17498-2921. The oscil- 
lations during big bursts at the known pulsar frequency show that these bursts were 
definitely from IGR J17498-2921. We find that at least several of the other bursts 
were also likely originated from IGR J 17498-2921. Spectral analysis reveals that the 
luminosity differences among various bursts are primarily due to differences in nor- 
malizations, and not temperatures, even when we consider the effects of colour factor. 
This shows burning on a fraction of the stellar surface for those small and medium 
bursts, which originated from IGR J17498-2921. The low values of the upper limits 
of burst oscillation amplitude for these bursts suggest a small angle between the spin 
axis and the magnetic axis. We find indications of the PRE nature of a medium burst, 
which likely originated from IGR J17498-2921. If true, then, to the best of our knowl- 
edge, this is the first time that two PRE bursts with a peak count rate ratio of as high 
as « 12 have been detected from the same source. 

Key words: methods: data analysis — stars: neutron — (stars:) pulsars: general — 
X-rays: binaries — X-rays: bursts — X-rays: individual (IGR J17498~2921) 



1 INTRODUCTION 

Periodic intensity pulsations at the neutron star spin fre- 
quency, thermonuclear X-ray bursts, and burst oscilla- 
tions (brightness variations close to the spin frequency 
during the bursts) are observed fro m some neut r on sta r 
low-mass X-ray binaries (LMXBs; iLamb et al.1 l|2009f ): 
IStrohmaver and BildsterTI jioOGI ')). Modelling of these fea- 
tures can be useful to measure the neutron star parameters, 
and to probe the strong gravity regime f e.g. . iBhattacharvva I 
l|2010t ):|Psaltis (200S))- However, in order to use these fea- 
tures as tools, one needs to understand them sufficiently 
well. Although, their basic properties are well understood, 
there are several outstanding questions. For example, why 
some of the neutron star LMXBs are accreting pulsars 
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(show i ng p e riodic pulsati o ns), w hile others are not (e.g., 
lOzel I (|2009l '): [Lamb et al. I (|2009h '). what creates burning re- 
gion asymmetry during decays of some thermonuclear bursts 
and g ives rise to burst decay oscillations (Bhattacharvy^ 
(|20ld 'l and references therein), what causes the plausible 
confinement of th e burning reg i ons a s indicated from tim- 
ing analysis (e.g., IWatts et al. I (|2008l )). etc. Accreting pul- 
sars showing burst oscillations are the ideal sources to ad- 
dress these questions. In this paper, we report the results 
of detailed timing and time-resolved spectral analyses of 
the thermonuclear bursts from the Rossi X-ray Timing Ex- 
plorer (RXTE) proportional counter array (PCA) field-of- 
view (FoV) of the recently discovered accreting pulsar IGR 
J17498-2921 with n e utron star spin frequ ency 400.99 Hz 
(|Papitto et al. II2OIII ). iLinares etal] (|201ll ) reported burst 
oscillations from this source. We confirm this feature with 
a much higher significance. We also show spectral indica- 
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tion of thermonuclear bursts from confined regions on tiie 
neutron star surface. 



2 DATA ANALYSIS AND RESULTS 

We analyze all the RXTE PCA data (32 obslds between 
August 13 and September 22; 146.496 ks exposure) of the 
2011 outburst of the accreting 401 Hz pulsar IGR J17498- 
2921. Twelve thermonuclear bursts are detected in the entire 
data (see §[3]for discussions on thermonuclear origin, and on 
the sources which could give rise to these bursts). We carry 
out timing analysis and time-resolved spectroscopy, mostly 
using Good-Xenon data files (~0.95 us time resolution), in 
order to study the nature of these bursts. 

The properties of all the bursts are given in Table [T] 
The bursts can be divided into three groups, based on the 
pre-burst level subtracted peak count rates (/peak: counts 
per second per proportional counter unit (PCU)): two big 
bursts with /peak ~ 3600 — 3700, three medium bursts with 
/peak ~ 240 — 300, and seven small bursts with /peak ~ 
45 — 110 (see Fig. [!}. However, this figure also shows that 
the shape and duration of bursts do not clearly change across 
the groups. 

Now we search for oscillations from all the bursts in the 
entire PCA energy range (using all active PCUs). We start 
with the August 16 big burst, for which burst oscillations 
were reported (see §[I|. The entire burst (above 5% of the 
peak count rate) is divided into 33 (= M) s egments of 1 
s each . The Leahy normaliz ed power spectra (|Leahv et al. I 
l|l983t ): lvan der Klis I (|l989l )) from all of them are averaged 
to obtain a power spectrum of 1 Hz resolution; and a range of 
±3 Hz from the known pulsar frequency (§[!]) is searched for 
a candidate peak. This is because the burst oscillation fre- 
quency do not shift f rom the neutron star spin f r equen cy by 
more than 3 Hz (e.g.. lStrohmaver and Bildsten I (|200y)). We 
find a candidate peak of « 4.17 Leahy power at 401 Hz. The 
probability of obtaining a power this high in a single trial 
from the e xpected no i se dist ribution (2A/ = 66 degrees 
of freedom: Iron der Klis I (| 19891 )) is « 5.88 x 10~^ Consid- 
ering a number of trials of 72 (= 6 x 12; six 1 Hz frequency 
bins are searched for each of 12 bursts), the significance of 
detection is « 4.1(t (estimated rms amplitude ~ 4.6±0.2%). 
This suggests that a further and stronger detection would 
be required for confirmation of burst oscillations from IGR 
J17498-2921. Therefore, we perform a similar timing anal- 
ysis for the August 20 big burst for M = 33. A candidate 
peak (« 7.04 Leahy power) appearing at 401 Hz has the 
single trial significance of 1 — 2.24 x 10^^". Considering a 
number of trials of 72 (as before), the significance of detec- 
tion is ~ 8.8a, which confirms burst oscillations from IGR 
J17498-2921 (see Fig.[2l note that this burst originated from 
IGR J17498-2921 (§[3])). If some of the 12 bursts were not 
originated from IGR J 17498-2921 (§[2|, then the significance 
of oscillations from both the big bursts would be hig her. The 
dyna mic power spectrum (Strohmayer and Markwardt I 
1 19991 ) of the August 20 big burst shows that the oscilla- 
tions appear intermittently during burst decay, and there is 
no significant frequency evolution (panel a of Fig. [3]). Panel 
b of Fig. [3] shows the rms amplitude evolution during the 
August 20 big burst. However, burst oscillations are not de- 
tected from any medium or small burst. The 3cr upper limits 



of rms amplitude are 4.4%, 5.0% and, 5.1%, for all medium 
and small bursts combined (Fig. [2]), for all medium bursts 
combined, and for all small bursts combined, respectively. 

Next we perform time-resolved spectroscopy of each 
burst after dividing them into time segments with suffi- 
cient counts. Note that each small burst has just one seg- 
ment to maintain enough statistics. From each segment, 
we create an energ y spectrum with dead time correction 
Ivan der Klis 1 1 19891 ) , and a background spectrum from the 
pre-burst emission |Bhattacharvva and Strohmaver I l2006l : 
lisallowav et al. IbOOa ). considering onlv the top lavers of all 
active PCUs. We fit each energy spectrum in 3 — 15 keV with 
a standa rd absorbed blackbody model (p habs*bbodyrad in 
XSPEC; [strohmaver a nd Bildstein (|2006l )) for a fixed neu- 
tral hydrog en column density A'h = 2.87 x 10"^^ cm~^ 
([Torres et al. 1120111 ). considering a systematic error of 1%. 
The model fits the spectra well, with reduced < 1-0 (de- 
grees of freedom = 28) for « 60% spectra, and between 1.0 
and 1.5 for almost all other spectra. The two big bursts show 
significant cooling in the decay portions, and the best-fit 
blackbody temperature and normalization (defined in the 
caption of Fig. [3| evolve in a correlated way. The tem- 
perature profile of each big burst shows two peaks, while 
a normalization peak coincides with the temperature min- 
imum between two peaks (see Fig. [3] for the August 20 
burst). This is a clear signat ure of a photosp h eric r adius 
exp ansion (PRE) burst (e.g.. iGallowav et al. I ||2008|): see 
also Linares et al.l l|201 j ): Chakrabortv and Bhattacharvva"] 
(|2011bl )). The medium bursts also show a cooling trend dur- 
ing decay (e.g.. Fig. 3)). This figure also shows a somewhat 
correlation between the best-fit blackbody temperature and 
normalization, which is indicative of the PRE nature of the 
August 19 medium burst, (see §[3] for a discussion). We find 
that the best-fit blackbody temperatures of all the bursts are 
consistent with each other (see Fig. [5] and §[3]). However, the 
corresponding best-fit normalizations are correlated with the 
burst fiuence (integrated energy) values, and both parame- 
ters increase roughly by two orders of magnitude from small 
bursts to big bursts (see Fig. [S]). Finally, in order to track 
the temperature evolution in a model independent way, we 
plot the burst colour (ratio of pre-burst level subtracted 
count rate above 6.14 keV to that below 6.14 keV) with 
time (Fig. [6|. This figure supports the finding that the tem- 
peratures of all the bursts are consistent with each other, 
and also shows a cooling trend during the decay of small 
bursts. A dip in colour (Fig. [6]) near the peak count rate 
supports a plausible PRE nature of the August 19 medium 
burst (see above). 



3 DISCUSSION AND CONCLUSIONS 

We discuss the implications of our results in this section. 
Let us first show that all the 12 bursts are thermonuclear 
bursts. The sharp rise and slow decay of intensity, accept- 
able fitting of burst spectra with a blackbody model, cool- 
ing during burst decay and detection of burst oscillations 
leave no doubt that the big bursts are of thermonuclear 
origin. The first three properties are also true for medium 
and small bursts. The very similar shape and duration of 
all the burst profiles argues that, if big bursts are ther- 
mormclear, then others from the same source are thermonu- 
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clear bursts too. Besides these are not repetitive bursts, 
and hence are not accretion-powered ty pe-II bursts (e.g., 
IChakrabortv and Bhattacharvval ((201l3))- These establish 
that all the 12 bursts are thermonuclear. 

The PCA FoV of IGR J17498-2921 contains five addi- 
tional known thermonuclear X-ray bursters: XTE J1747-274 
(or, IGR J17473-2721), SAX J1750.8-2900, SAX J1747.0- 
2853, lA 1742-289 (or, AX J1745.6-2901) and SLX 1744- 
299/300. This brings the question whether some or all of 
the 12 detected bursts (§ ^ were not from IGR J17498- 
2921. The two big bursts showed oscillations at the known 
pulsar frequency of 401 Hz, and hence they were definitely 
from IGR J17498-2921. But which sources gave rise to the 
medium and small bursts? The first four of the five above- 
mentioned additional bursters are transients, and none of 
them were in outburst during the outburst of IGR J17498- 
2921 (otherwise, t he scanning programs of satellites, such as 
INTEGRAL fe.g.. lGibaud et al. I (j2oTll ')') would detect such 
an outburst). Since a thermonuclear burst from the quies- 
cent p hase of a transient is very rare (e.g.. iKuulkers et al. I 
(|2009l )). it is unlikely that the medium and small bursts 
originated from one of the four known transient bursters, 
or any other heretofore unknown transient burster in quies- 
cence in the FoV. Here we consider that the medium and 
small bursts have originated from IGR J17498-2921 and/or 
from the persistent burster SLX 1744-299/300 though there 
still remains a small chance that a low-intensity outburst 
(missed by X-ray satellites) of a transient at the edge of the 
FOV could give rise one or more of these bursts. Three bursts 
were previously observ ed with PCA from SLX 1744-299/300 
l|Gallowav et al. 1 120081 '). We find that the shape and dura- 
tion of these bursts were similar to those of our 12 bursts. 
The peak count rates (when corrected for the off-axis po- 
sition of the source) of all these SLX 1744-299/300 bursts 
are less than those of our medium bursts, but greater than 
those of our small bursts. Therefore, the shape, duration 
and peak count rate do not reveal whether our small bursts 
and medium bursts were from SLX 1744-299/300 or from 
IGR J 17498-2921. The Table 3 of iGallowav et al. I (|2008l ') 
shows that the occurrence of SLX 1744-299/300 bursts per 
hour was 0.037 for 290 ks of PCA observation. If all our 
medium and small bursts had originated from SLX 1744- 
299/300, then the burst rate of this source would be 0.246 
hour~^, which is about seven times larger than the known 
burst rate (0.037 hour~^). Even if only the small bursts 
had come from SLX 1744-299/300, its burst rate would be 
about five times larger than the known rate. These suggest 
that not all the medium and small bursts may have origi- 
nated from SLX 1744-299/300, which imphes that at least 
some of these bursts were from IGR J17498-2921. In or- 
der to know whether this is indeed the case, we attempt 
to find out if the angular location of a burst is consistent 
with that of IGR J17498-2921. For this, we e xploit the fact 
that the 5 PCUs are not perfectly aligned (|jahoda et al. I 
l2006l ). and hence the ratio of observed count rates in a pair 
of PC Us depends upon the p osition of the source within the 
FoV l|Gallowav et al. Il2008f ). Two PCUs were on for each 
of the bursts, except for the August 16 big burst (1 PCU). 
Therefore, for each of 11 bursts, we compute the ratio {Ri, 
with an error) of observed total counts (preburst level sub- 
tracted and deadtime corrected) in the pair of PCUs. For 
the PCA observations of IGR J17498-2921, this source was 



almost (within a few arcsec) at the centre of the FoV. There- 
fore, if an Ri value is consistent with the expected value for 
the centre of PCA FoV, then the corresponding burst likely 
originated from IGR J17498-2921. In order to find out these 
expected values, we consider several bursts in PCA observa- 
tions from two other sources (4U 1636-536 and 4U 1608-52), 
which were at the centre of PCA FoV, and for each of which 
there is no known burster in the PCA FoV. For these sources 
and for each pair of PCUs, we compute a mean burst count 
ratio {R2, with an error), in the same way we compute Ri. 
These _R2's are the expected ratio values for the centre of 
PCA FoV. We compare an Ri value with the R2 value for 
same PCU pair. For each of 11 bursts, Table [1] gives the 
difference between Ri and R2 in the unit of the estimated 
error in this difference. A smaller value of this difference im- 
plies a higher possibility for a burst to be originated from 
IGR J17498-2921. For example, this difference is less than 
one for the brightest medium burst (August 19; for which 
a plausible PRE nature is seen (§[2])) and two small bursts, 
implying that these bursts were likely originated from IGR 
J 17498-2921. 

Based on the above discussion and Table[T] in the rest of 
this section we assume that at least one medium burst and 
at least one small burst were originated from IGR J17498- 
2921. We now ask the question what made some bursts 
from IGR J17498-2921 so energetic, while one or more other 
bursts were so weak. Was it because of burning of different 
chemical compositions in differe nt accretion rate regimes 
(|Strohmaver and Bildstenl (|2006l ) and references therein)? 
We find that this is very improbable, because (1) the shape 
and duration of all the burst profiles are similar, indicat- 
ing the burning of similar compositions; and (2) the burst 
fluence is not correlated with the persistent intensity. For 
example, the persistent intensities before one of the August 
19 small bursts and before the August 29 small burst are 
respectively slightly higher and slightly lower than the per- 
sistent intensity before the August 20 big burst. Hence we 
try to find if the increase of the fiuence from small bursts 
to big bursts for IGR J17498-2921 is due to the blackbody 
temperature increase or due to the burning area (oc black- 
body normalization) increase. Figs. [3] (c), [4] (6), [5] (a) and[6] 
show that the best-fit temperatures of one burst are consis- 
tent with those of others. In Fig. [5] (a), the temperatures of 
small bursts are somewhat systematically less than those of 
medium and big bursts; but this is because the former tem- 
peratures are average values during burst profiles, while the 
latter temperatures are maximum values during burst pro- 
files. Note that such a maximum temperature value for a big 
burst is for the second temperature peak, which corresponds 
to the touch-down point (i.e., settling of the photosphere on 
the stellar surface after expansion). Unlike the temperature, 
the normalization, which is proportional to the burning area 
(see caption of Fig. [3]), increases significantly with fluence 
(Fig.[S](6); §[2]). This strongly suggests that the burning area 
increases from small bursts to big bursts for IGR J17498- 
2921, while the other parameters, including the temperature 
and composition, remain roughly unchanged. If this is true, 
then each burst from IGR J17498-2921 should be similarly 
strong within its confinement (burning area), and hence, 
even the small and medium bursts from this source could 
have local photosphere expansion. However, for a smaller 
burst, a signature of such expansion (discussed in § [2| is 
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washed away due to the large time bins required to gather 
sufficient counts for spectral analysis. Nevertheless, such a 
signature appears for the brightest medium burst (Figs. 3] 
[b, c); §[11), which likely to be originated from IGR J 17498- 
2921 (see Table [T] burst 3). Although, a temperature peak 
before the peak intensity is not there, the following proper- 
ties indicate the photosphere expansion: (1) normalization 
evolution is very similar to that of the big bursts (compare 
Fig. (d) with Fig. 2] (c)); (2) intensity peak corresponds 
to the normalization peak and a low temperature; and (3) 
a temperature peak (corresponding to a low normalization) 
appears after the intensity peak. Moreover, an observed dip 
in colour near the burst intensity peak (Fig. [6](c)) suggests 
a similar dip in temperature (undetected plausibly due to 
large time bin as mentioned above), and hence a plausible 
temperature peak before the peak intensity. However, a con- 
clusion about the PRE nature of this burst has to be made 
cautiously, because (1) a temperature peak before the peak 
intensity is not significantly detected, and (2) the normaliza- 
tion value following the maximum is not lower with at lea st 
4a significance (as considered by iGallowav et al. I l|2008l )). 
If this medium burst is really a PRE burst, then, to the 
best of our knowledge, this is the first time that two PRE 
bursts with a peak count rate ratio of as h igh as ~ 12 have 
been d etected from the same source (e.g., IGallowav et al. I 
l|2006bD ). This makes the standard method of source dis- 
tance measurement using PRE bursts fe.g.. iKuulkers et al. I 
l|2003l )) somewhat less reliable, at least for pulsars. All the 
bursts from IGR J17498-2921 were likely helium-rich, be- 
cause of there high temperatures, somewhat short durations 
and the PRE nature (for at least the tw o big bursts; see, for 
example, IGallowav et al. ] (|2006al . l2008l )). 

The blackbody normalization values for small and 
medium bursts (Fig. [5] (6)) imply burning areas much 
smaller than any realistic neutron star surface area (e.g., for 
a stellar radius of 8 — 20 km) , implying burning in confined 
regions. Note that these normalization values are estimated 
assuming a source distance of 10 kpc (see caption of Fig. [3]). 
If the distan ce is less (e.g. , 7.6 kp c for IGR J17498-2921 as 
reported bv lLinares et al.l (|201ll )l. the intrinsic normaliza- 
tion values will be smaller (implying even smaller burning 
areas) by a factor same for all the bursts from a given source. 
A correction due to t he surface gravitational redshift 1 + z 
l|Sztaino et al. iri985h will make the intrinsic normalization 
values further smaller by another factor (1 -I- z)^ , which is 
same for all the bursts from a given source. 

The absorption and scattering in the neutron star 
atmosphere makes the observed temperature higher (rel- 
ative to the intrinsic temperature) by a factor /, and 
the observed normalization lower (relative to the in- 
trinsic norma lization) by a factor f^, where / is the 
colour factor ([Sztaino et al. Ill985l: iMaiczvna et ahl l2005l : 



ISuleimanov et ahT 201ll :" Bhattacharvva et al. I 2010l ). Does 
this mean that our conclusion about burning in confined 
regions for IGR J17498-2921 is not robust? In order to 
find out, let us examine, if the burning areas of the big 
bursts cover the entire neutron star surface, then, whether 
the burning area of a medium burst from IGR J 17498-2921 
can also cover the entire surface. How could this happen? 
Suppose, the intrinsic temperature of the medium burst is 
smaller than that of big bursts, and / increases with the de- 
crease of temperature to make the observed temperature of 



the medium burst and big bursts similar. Then the ratio of 
a big burst intrinsic normalization to the medium burst in- 
trinsic normalization will be lower (by a factor, say, g) than 
the ratio of a big burst observed normalization to a medium 
burst observed normalization (the observed normalizations 
are given in Fig. [S] (6)). If g ~ 10, then the burning ar- 
eas of the medium bursts could be similar to those of big 
bursts (Fig. [5] (6)). But, considering the relevant extreme 
/ valu es (1.64, 1.22) from the tables of IMaiczvna et al.l 
hoO^ ). g < 3.3(= [1.64/1.22]"). Therefore, even if we con- 
sider the maximum possible value (i.e., the entire surface) 
for a big burst burning area, a medium burst (from IGR 
J17498-2921) burning area cannot be more than ~ 35% of 
the stellar surface, even after considering the effects of colour 
factor. With the same arguments, a small burst (from IGR 
J17498-2921) burning area will be even smaller. This shows 
that our conclusion about burning in confined regions is ro- 
bust. Finally, we note that the observed burning areas of 
some of the bursts may be small due to obscuration (for 
example, if the observer is close to one spinning pole, and 
burning regions are close to the other); but even for this, 
the burning has to happen in confined regions. Even if all 
the medium and small bursts were from SLX 1744-299/300, 
the difference in normalization values between these two sets 
of bursts would indicate confined burning for small bursts 
for this source. However, we note that the confined burning 
could be more easily explained for IGR J17498-2921 (than 
for SLX 1744-299/300), which is a known pulsar, and hence 
likely to have a higher neutron star magnetic field (see the 
next paragraph). 

Our results support a prediction of iBildsten I (|l993f ). 
that the weaker bursts are from small fractions of the neu- 
tron stars. The low values of the upper limits of burst os- 
cillation fractional rms amplitude for small and medium 
bursts (see § [2| suggest that the burning regions for IGR 
J17498-2921 were close to a spi nning pole (assuming a hot 
spot model; iLamb et al. I (|2009l )). Hence, a reasonable as- 
sumption that the thermonuclear burning happens close to 
the magnetic pole(s) (as the accreted matter is channeled 
to these poles for a pulsar) implies that the spin axis and 
the magnetic axis of this pulsar are close to each other. 
This supports iLamb et al. ] ((200i)'s model of accreting mil- 
lisecond X-ray pulsars. What gave rise to oscillations dur- 
ing the decays of big bursts? If the burning covered the 
entire neutron star surface, then the oscillations could be 
due to a complex asymmetric brightness pattern, caused by 
Rossby waves (or r-modes) in the surface lay er s, shear insta - 
bilities, etc . (e.g iLee fc Strohmaver I (120051 ): iHevl I (|2005h : 
ICumming I (|2005l ): |Piro and Bildsteni (|2006D '). Finally, what 
confined the burning region, and why did the burning area 
change from burst to burst? Detailed theoretical and nu- 
merical studies are required to answer these questions. The 
confined burning could be due to the m agnetic field near the 
magnetic pole (but see IWatts et al. I (j2008)), or, magnetic 
field locally enhanced by th e convection durin g thermonu- 
clear flame spreading (e.g., ISpitkovskv et al.l ((2002 ') )■ or, 
mag netic structure (e.g., belt) o n the neutron star surfaces 
re.g.. |Pavne and Melatosl (|2006l ')'). or, higher local tempera- 
ture (sufficient for ignition) near the magnetic poles, while 
ignition condition is not reached at rest of the stellar sur- 
face (Watts et al. llioO^ '). or, something else. Whatever be 
the reason, the indication of burning in confined regions 
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reported in this paper may have significant impact in un- 
derstanding the accreting pulsars, burst oscillations, burst 
ignition conditions, surface magnetic field structure and its 
interaction with flame spreading, etc. (see also §[!]). 
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Figure 1. Example of each of three types of bursts, big (left 
panel; August 20, 2011), medium (middle panel; August 19, 
2011) and small (right panel; August 29, 2011), from the IGR 
J17498-2921 observations. Each RXTE PCA light curve has 
1 s time binning. Note that the peak count rate of the big 
burst is about an order of magnitude larger than that of the 
medium burst, and the latter one is a few times larger than 
the peak count rate of the small burst. In the insets, nor- 
malized profile of the August 20 big burst is compared with 
that of three combined medium bursts, and that of seven 
combined small bursts, after aligning the peaks. These in- 
sets (with 3 s time bins) show, despite a large change of peak 
count rate from one burst to another, the shape/duration of 
the big, medium and small bursts are similar to each other 
(see §[2ll. 
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Figure 2. Leahy normalized power spectra from the RXTE 
PCA data of IGR J17498-2921. Panel a: power spectra from 
33 time segments (Is each) during the big PRE burst (Au- 
gust 20, 2011) are averaged without any further frequency 
rebinning. The strong peak at 401 Hz shows significant burst 
oscillations. Panel b: same as panel a, but the power spectra 
from all small and medium bursts are merged together. The 
power spectrum, which is an average of 331 power spectra, 
does not show burst oscillations (see § [2]| . 
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Figure 3. Properties of the big burst (August 20, 2011) 
with oscillations from IGR J17498-2921. Panel a: RXTE 
PCA light curve (soHd curve) and dynamic power spectrum 
(contours). The latter shows 30% to 99% contours of the 
maximum power 37.8, using 2 s overlapping time bins with 
0.25 s shift between two adjacent bins. These contours sug- 
gest the presence of an 401 Hz signal with no significant 
frequency evolution. The dotted horizontal line shows the 
pulsar frequency (see § [2}. Panel b: fractional rms ampli- 
tudes in 1 s bins during the burst, la error bars are given 
for 3cr significant four amplitudes (6.24±0.39%, 7.14±0.38%, 
11.53±0.53% and 11.23±0.66%), and 3a upper limits are 
given for the rest (see § [Jjl. Panel c: the evolution of the 
best-fit (in 3 — 15 keV) blackbody temperature (with 90% 
errors) of the burst (see § O. The rise of the temperature 
is not seen because of the somewhat large time bin. Panel 
d: the evolution of the best-fit (in 3 — 15 keV) blackbody 
normalization (with 90% errors) of the burst. This normal- 
ization is proportional to the burning area, and is defined 
as ^|„/^ioi where -R^m is the neutron star radius in km 
when the entire surface emits, and Dio is the distance to the 
source in the unit of 10 kpc. The specific correlation between 
the temperature and the normalization shows that this is a 
PRE burst (see §[2]l. 
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Figure 4. Parameter evolution for a medium burst (August 
19, 2011) from an IGR J17498-2921 observation. Panel a: 
count rate evolution with 1 s time bins. Panel b: the evolution 
of the best-fit (in 3 — 15 keV) blackbody temperature (with 
90% errors) of the burst. Panel c: the evolution of the best- 
fit (in 3 — 15 keV) blackbody normalization (defined in the 
caption of Fig. [S] with 90% errors) of the burst. The specific 
correlation between the temperature and the normalization 
(especially when compared with panels c and d of Fig. [Sjl 
indicates photospheric radius expansion, although the first 
temperature peak is not visible (see § [2] and [3]l . 



Table 1. Properties of bursts from the 2011 RXTE PGA observations of IGR J17498-2921. 



Serial 
no. 



Observation 
start time^ 



Time in 
MJD2 



Peak count 
rate^ 
(count /s) 



Temperature® 
(keV) 



Normalization" 



\Rl -B.2\/a, 



7 



{R1-R2) 



(s) 


6.07 




114.351}*:^° 




4.13 


r, 07 + 0.49 

^••^'-0.37 


0.87«:I1 


1.37 


5.20 


77+0.35 
'-0.28 




0.65 


4.17 


-I 7-1+0.35 
^- ' ^-0.28 


^■^'^-1.29 


1.03 


4.31 


1 QQ+0.34 


1.0»_Q 7g 


1.01 


4.12 


r. qo + 0.26 

z.oo_g 22 




6.89 


7.05 


r, cn + o.oe 
^•'^^-o.oe 


111.20+11 


1.72 


4.17 


r)Q + 0.28 

^•^"-0.24 


8-85_2;93 


1.96 


3.30 


z.o»_Q 29 


n 70+0.39 


0.88 


6.20 




1 40+0-95 
^■^^-0.60 


1.15 


6.06 


or\+0.35 
^•^"-0.28 


^- '°-0.32 


0.26 


3.13 


, 7c;+0.51 
^•"^-0.38 




2.93 



1 


2011 


08- 


16T15 


19 


28 


55789 


64 


2 


2011 


08- 


19T12 


11 


28 


55792 


53 


3 


2011 


08- 


19T13 


45 


20 


55792 


58 


4 


2011 


08- 


19T15 


19 


28 


55792 


67 


5 


2011 


08- 


20T10 


08 


32 


55793 


43 


6 


2011 


08- 


20T10 


08 


32 


55793 


44 


7 


2011 


08- 


20T13 


16 


32 


55793 


59 


8 


2011 


08- 


21T12 


46 


24 


55794 


57 


9 


2011 


08- 


29T08 


47 


28 


55802 


39 


10 


2011 


09- 


03T07 


50 


24 


55807 


34 


11 


2011 


09- 


08T05 


49 


20 


55812 


25 


12 


2011 


09- 


22T09 


25 


20 


55826 


41 



3648.72±61.32 
109.04±14.90 
291.08±20.05 
58.08±12.96 
76.69±13.89 
242.85±18.98 

3634.12±61.20 
273.42±19.47 
84.51±13.31 
60.40±11.96 
45.28±11.06 
83.99±11.88 



The start time of the event data file containing the burst. 
The burst peak occurence time in MJD. 

The preburst level subtracted burst peak count rate per PGU (PGU2), Itr error is given. 
The exponential decay time for each burst. 

The best-fit blackbody temperature in 3 — 15 keV, 90% error is given (see Fig.fSj. 

The best-fit blackbody normalization in 3 — 15 keV, 90% error is given (see Fig. [S] defined in the caption of Fig. [3t. 

R\: ratio of observed total counts of a burst in the pair of PGUs; -R2: expected value of Ri, if the burst were from IGR J17498-2921 
(see §[3ll. c!(Fi-^^—ii2) estimated Ict error in _Ri — R2. Hence a lower value of |_Ri — R2\/ implies a higher possibility for a 

burst to be originated from IGR J17498-2921. For the burst 1, only one PCU was on. 
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Figure 5. Best-fit blackbody temperature (panel a) and 
normaUzation (panel 6; defined in the caption of Fig. [Sj in 
3—15 keV as functions of tfie burst fluence for all RXTE 
PC A bursts from IGR J 17498-2921 observations. Note that 
for each small burst, only one spectrum for the entire burst 
duration is fitted, and the corresponding temperature and 
normalization are plotted. For each medium and big burst, 
the maximum temperature reached during the burst and its 
corresponding normalization are plotted. Cross signs: big 
bursts; square signs: medium bursts; triangle signs: small 
bursts. 90% error bars are given. Panel a: the dotted hor- 
izontal line goes through the big burst (August 20, 2011) 
temperature. Panel b: the dotted horizontal lines from top 
to bottom go through the big burst (August 20, 2011) nor- 
malization, mean normalization of medium bursts, and mean 
normalization of small bursts respectively. This figure shows 
that, while the temperatures of all the bursts are similar to 
each other, the normalization (and hence the burning area) 
is correlated with the burst fluence, and it increases from 
small bursts to big bursts by about two orders of magnitude 
(see §[2]l. 
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Figure 6. The colour (cross sign witli Icr error; ratio of 
the persistent-subtracted count rate above 6.14 keV to that 
below 6.14 kev) and count rate (solid curve) evolution dur- 
ing the bursts from the RXTE PCA observations of IGR 
J17498-2921. Panel a: big burst (August 20, 2011). Panel 
b: three medium bursts averaged. Panel c: the brightest 
medium burst (August 19). Panel d: seven small bursts aver- 
aged. Burst peaks are made aligned for averaging. This figure 
shows that (1) the maximum values of colour for big, medium 
and small bursts are similar to each other; (2) there is a de- 
creasing trend of colour during the decays of big, medium 
and small bursts; and (3) a dip in colour near the peak count 
rate of the brightest medium burst (August 19) is consistent 
with a plausible PRE nature (see §[2j. However, we note that 
such a dip is not seen in the average profiles (panels b and 
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